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ABSTRACT 

In a frequency domain in which a material thickness is smaller than the skin 
depth, a formalism based on the sheet impedance concept has been developed 
and introduced in the FD-TD code ALICE. The predictive capabilities of the 3D 
code have been evaluated by comparison to analytical and experimental data. 

1. INTRODUCTION 

For large structures such as aircraft, the spatial resolution of 3D computer 
codes is insufficient to model thin surfaces such as lossy skins. There is 
therefore a need to accurately include their effects on the electromagnetic 
coupling paths and so, to define an approach whereby lossy surfaces can be 
treated in the time-domain. This can be achieved by defining a relationship 
between the surface electric and magnetic fields such that the details of the 
skin do not have to be directly modelled. In a frequency domain in which a 
material thickness is smaller than the skin depth, a formalism based on the 
sheet impedance concept has been developed and introduced in the FD-TD code 
ALICE [1]. The predictive capabilities of the 3D code have been evaluated by 
comparison to analytical results obtained by K.F. Casey [2] on a loaded 
aperture and by comparison to experimental data obtained during an in-flight 
experiment conducted in France to study lightning interaction with a vehicle 
[3]- 


2. LOW FREQUENCY ELECTROMAGNETIC PENETRATION OF LOADED APERTURES 


In a frequency domain in which a material thickness is smaller than the skin 
depth, it can be shown [4] that a lossy material may be characterized by its 
sheet impedance defined as: 
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where J is the surface 

8 

field component. 


current density and E is the tangential electric 
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For good conductors in 
compared with conduction 
equation (1) reduces to: 


which displacement currents are negligibly small 
currents for the frequency of interest (cr >> £ou) , 


1 



where cr and d are the material conductivity and thickness respectively. 
The Magnetic Shielding Effectiveness (MSE) of a structure is defined as: 


MSE 



(3) 


where H # is the field at a point chosen in the absence of the shield and H is 
the field at the same point with the shield in place. 


K.F. Casey [2] has shown analytically that the MSE of a circular aperture (of 
radius a) loaded by a panel characterized by its sheet impedance Z is given, 
in the frequency domain, by: 


with: 
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R: the net resistance of the joint (length L) located around the panel. 


F^S 


F t is obtained by solving the linear system F^ 

3n F 

3 = j and K is a combination of m and n index. 
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the metallic plane is perfect 

(7) 

This expression corresponds to a first order filter function with a cut-off 
frequency F - F for 3 dB attenuation. 

C CO 

3. FD-TD MODELLING 


Vhen the junction between the panel and 
(R =0), equation (4) reduces to: 


MSE = 


1 + 


F 

j — 
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If we consider a composite sheet in the plane OXY of a unit cell (figure 1), 
it can be shown [5] that the current density J may be expressed as: 


1 

J : E (8) 

Z Oz 
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With this expression for J , equation (2) in reference [1] can be solved in 
finite difference form without taking into account the panel thickness. 
Assuming that the contact resistance around a panel corresponds to a joint of 
width W, we can define a surface impedance of this joint as: 


RL 

v r 

The expression (6) for F c then becomes: 
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In the FD-TD code, W will be the mesh size. 

The joint surface impedance Z^ can be related to a conductivity and a 
thickness dj so that (o^d^)* 1 . 

4. COMPARISON TO ANALYTICAL RESULTS 

The structure used to study field penetration through a loaded aperture is a 
perfectly conducting cavity having a square aperture on the upper side 
(figure 2). The aperture dimensions are small compared to the cavity size. 
The box is connected to two thin wires of infinite extend in space. A current 
generator is located on one of them at 5.3 m from the cavity. The waveform of 
the injected current is of trapezoidal shape with a rise time of 2 ns and a 
peak value of 1000 A. In a first step, the joint around the panel will be 
considered as perfect (R = 0) . 
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Figure 3a shows the waveform of the magnetic field component H # computed at a 
point P (figure 2) in the case of a free aperture. 

The point P is located inside the cavity, 5 cm below the aperture. As 

indicated by the FFT of H o shown in figure 3b, the high frequency content is 

mainly due to a X/4 resonance of wire 1 (between the current generator and 

the box - F = 15 MHz) and to a \/2 resonance of the cavity (F = 45 MHz) . 

Figures 4a and 5a show the temporal waveforms of the magnetic field component 
H computed at the point P when the aperture is loaded with a composite panel 
having a sheet impedance of 10 and 0.1 0 respectively. A low frequency 
content appears as a consequence of the panel filtering effect. 

Figures 4b and 5b show the FFT ratios for the magnetic field components H and 
H o corresponding to the previous results. The agreement with the theory 
presented in S 2 is excellent. The code has predicted: 

- The same low frequency content for H o and H. 

- The cut-off frequency corresponding to a 3 dB attenuation (7), (5). An 

equivalent radius "a" of the square aperture can be calculated by 
comparison to a circular aperture having the same area. Using (5) with 

a : 0,3 i, we obtain F = 1 MHz for 1 - 1 ft and F = 0.1 KHz for Z = 0.1 ft. 

CSC • 

- The decrease of the spectrum with a slope of 20 dB/decade for F >> F fi as 
indicated by (7) . 

We consider, now, that there is a joint all around the panel. Two values 

Z i - 2.5 ft and Zj = 22.5 ft have been considered. These values have been chosen 

to satisfy F /F =2 and F /F = 10. With the contact length L = 2 tta, the 
corresponding contact resistances are R = 0.16 ft and R = 1.43 ft respectively. 
From (5), the theoretical cut-off frequency for a loading material having a 
sheet impedance Z # = 1 ft is F c# = 1 MHz. 

Figures 6 and 7 show the variation of the MSE given by (4) and (7) and also 
the variation of the MSE obtained with the FD-TD code using W = 0.12 a: H # 
and H are the magnetic field components computed on the axis of the aperture, 

15 cm inside the cavity, in the case of a free aperture and in the case of a 

loaded aperture with a joint, respectively. For R = 0.16 ft, the approximate 
expression (7) is close to the exact value given by (4) and the FD-TD result 
is also in very good agreement with these results. For R = 1.43 ft, the 

cut-off frequency (for 3 dB attenuation) is 6 MHz with the exact expression 
(4) and 10 MHz with the approximate expression (7). The FD-TD result gives a 
cut-off frequency of 7.8 MHz between the two previous values. 

5. IN-FLIGHT EXPERIMENT MODELLING 

In reference [1] when the current pulse simulating the negative leader 

mechanism is injected onto the aircraft by means the 3D computer code all the 
calculated fields are in good agreement with experimental results except the 
magnetic field computed in the fuselage, behind the carbon composite door. 
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The signal obtained (figure 8b) exhibits a first order low-pass filter 
response having a cut-off frequency of 9 kHz whereas the measured signal 
(figure 8a) exhibits a response with a 50 kHz cut-off frequency- It has been 
shown that the only possible path for this high frequency content is the 
joint located all around the door and used to pressure the aircraft: a joint 
of poor quality increases the penetration of electromagnetic energy and so 
increases the internal frequency content (see § 2) . This joint resistance 
being unknown our purpose here is to calculate its value and then to 
introduce this value in the 3D code and compare the results to the 
experimental data. 

When penetrating inside the fuselage the external magnetic field is: 

- Attenuated by a geometrical effect. This attenuation (40 dB) has been 
evaluated using the 3D code and assuming the door is unloaded (figure 9) . 

- Filtered by the carbon composite material of the door. This filtering 
effect if given by the relations (4) and (6) . 

Knowing the values of F c , a and Z ( , we can deduce from (6) the resistance R 
of the joint. In the case of the Transall experiment, we obtained Z f = 30 mQ, 
a = 0.68 m and F c = 50 kHz, so R = 12 aft. 

These values of Z and R have been introduced in the code "ALICE" and the 
result is shown in figure 10. This response is in very good agreement with 
the in-flight measurement. 


6. CONCLUSION 

The goal of this paper is to present a formalism based on the sheet impedance 
concept which has been introduced in a 3D FD-TD code in order to take into 
account electromagnetic coupling through lossy materials and through 
resistive joints. The predictive capabilities of the code have been evaluated 
by comparison to analytical results obtained by K.F. Casey. The numerical 
tool has been used to model the penetration of electromagnetic fields through 
a carbon composite door located on the fuselage of an aircraft struck by 
lightning during an in-flight experiment. A very good agreement with 
in-flight measurements can be obtained if one takes into account the coupling 
through the joint surrounding the door. 
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Figure 1 : Elementary cell with a 

composite material (e.o.d). 


Cavity 


Aperture (0.54 x 0.54 ■) 



Figure 2 : Structure used to calculate field penetration 
through a loaded aperture. 


43-6 




A/m 


0. 0-2 0.4 0.6 0.6 1.0 1.2 1 4 1.6 18 2.0 2.2 2 4 2.6 2.S 3.0 

T(1E-6 s) 


111! 




■Oil 


F(Mhz) 45 


Figure 3a : Waveform of the magnetic component Figure 3b : FFT of H # component, 
at point P. {free aperture) 
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Figure 4a : Waveform of the magnetic component Figure 4b : B and H^FFT ratio, 
at point P. 

{loaded apertur.e Z § = 1 fl) 
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Figure 5a : Waveform of the magnetic component Figure 5b : H and H # FFT ratio 
at point P. 

(loaded aperture Z.= 0.1 fl) 










Figure 6 


Magnetic field attenuation 
for a composite panel 
(Z # = 1 £1) with a net contact 
resistance R = 0.16 fi. 
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Figure 8 : Internal magnetic field. 

a) Measured during in-flight experiment. 

b) Calculated without taking into account 
the joint resistance. 
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Figure 9 : Calculated magnetic fields. 

a) On the external skin of the fuselage 
near the composite door. 

b) Inside the fuselage assuming the door 
unloaded. 



Figure 10 : Magnetic field amplitude behind the 
composite door of the aircraft. 
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